Finite-temperature lineshapes in gapped quantum spin chains 
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We consider the finite-temperature dynamical structure factor (DSF) of gapped quantum spin 
chains such as the spin one Heisenberg model and the disordered transverse field Ising model. At 
zero temperature the DSF in these models is dominated by a delta-function line arising from the 
coherent propagation of single particle modes. Using methods of integrable quantum field theory 
we determine the evolution of the lineshape at low temperatures. We show that the line shape is in 
general asymmetric in energy and we discuss the relevance of our results for the analysis of inelastic 
neutron scattering experiments on gapped spin chain systems such as CSN1CI3 and YBaNiOs. 



Quasi one-dimensional spin chains are materials where 
quantum fluctuations give rise to striking strongly corre- 
lated phenomena. An exemplar of such behavior is the 
distinction, first identified by Haldane 1 - over 20 years ago, 
between integer and half-integer isotropic spin chains. 
The former are generically gapped while the latter are 
generically gapless. This effect is topological in origin 
and arises from the presence of a quantized Berry's phase 
in the effective model describing the chains. Contempo- 
rary examples of strong correlations in spin chains revolve 
around the role the multi-excitation continuum play in 
their physics. This continuum has been studied both 
experimentally-'- 3 -^ and theoretically^ and is understood 
to be the origin of a process known as spectrum ter- 
mination, where coherent excitations cross into a multi- 
excitation continuum and then experience rapid decay^ 

In order to probe the dynamical behavior of spin 
chains, inelastic neutron scattering is the premier tool. 
In particular, using inelastic neutron scattering it is pos- 
sible to determine with impressive accuracy the spectrum 
of spin excitations together with their lifetimes^ The 
theoretical counterpart of such measurements is a com- 
putation of the dynamical structure factor (DSF). For 
theoretical models that admit an integrable continuum 
field theoretic description, the computation of the zero 
temperature DSF is possible with impressive accuracy*^ 
However at finite temperatures, while important progress 
has been made /^ 11 ' 12 a general theoretical framework 
has yet to be settled upon. It is one aim of this letter to 
outline a promising approach. 

We do so against an experimental mise en scene of par- 
ticular relevance. In a system that supports a coherent, 
gapped, magnetic single-particle excitation at T = 0, the 
question arises of how the corresponding delta- function in 
the DSF broadens at non zero-temperatures i 13 ' 14 i 15 ' 16 i 17 
A partial answer to this question has been given by 
Sachdev and collaborators! 13 > 14 they demonstrated that 
at temperatures far below the gap the broadening in the 
immediate vicinity of the T — gap is Lorentzian in form. 
In the present work, using our approach to computing 
finite temperature DSFs, we determine the entire line- 
shape. As our central finding, we demonstrate that the 
lineshape is always asymmetric in energy, a feature that 
becomes more pronounced as the temperature increases. 



While we focus upon the lineshape in gapped quantum 
spin chains, we stress that this approach is applicable to 
the computation of general response functions in generic 
inte grab le continuum models, such as those considered in 
Ref.JlllIl). 

General Theoretical Framework: The systems we 
study all have representations as general Heisenberg 
models: 

H = 2^ ^-J-S_i_i • Sj_ i+1 + J z S Z iS Z i + \ + H • S. (1) 



Here Si = (Sxi, S z i) is a quantum spin (either integer or 
half-integer) at chain site i. We allow both the spin-chain 
to have anisotropic couplings ( Jj_, J z ) and for a magnetic 
field, H, to be potentially present. We are interested in 
computing the DSF 
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To compute this quantity, we expand C(t,x) = 
(S(r, x)S(O)) in a basis, {\l)}, of exact eigenstates of H, 

C(t,x) = ^^2e-^(l\S z (T,x)\m)(m\S z (0)\l), (3) 

l,m 

where E\ is the energy of eigenstate, and Z — 



-0E, 



is the partition function of the theory. By 



virtue of the gap, A, in this system, the Fourier transform 
of C(r, x) has a well defined low temperature expansion. 

This representation of the DSF finds its virtue when 
we employ a continuum, integrable reduction of the (lat- 
tice) model in Eqn. ([T]). In such cases the matrix ele- 
ments {l\S z (0)\m} can readily be computed exactly. At 
T = this permits the exact computation of Imx(^,q) 
at energies, u>, in the vicinity of the gap through the 
computation of a small number of matrix elements. At 
finite temperatures, this approach, for the problem at 
hand, breaks down in two fashions: i) the needed ma- 
trix elements (as well as Z) become highly singular ob- 
jects; and ii) to obtain the finite temperature broadening 
of the coherent mode, an infinite number of matrix ele- 
ments are needed. We solve these problems in a two step 
fashion. The singularities of the matrix elements are inti- 
mately associated with treating the spin chain as infinite 
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in length. While it is possible in certain circumstances 
to deal with these singularities directly ^ 10 ' 11 ? 12 ! 19 to cir- 
cumvent this first difficulty, we instead work with chains 
of large but finite length, R. The infinities in the matrix 
elements are then reduced to terms merely proportional 
to R which are cleanly cancelled by similar terms in the 
partition function. As part of this, we will exploit the 
fact that the matrix elements, even at finite R 7 are com- 
putable up to exponentially small corrections. To handle 
the second difficult, we recognize that the infinite sub- 
set of needed matrix elements from the sum in Eqn. Q 
are organized according to a Dyson's equation. This al- 
lows us to characterize the subset by resumming a finite 
number of matrix elements. We now consider how this 
approach works in practice in two experimentally rele- 
vant cases, the transverse field Ising model and the spin-1 
chain as represented by the 0(3) non-linear sigma model. 
Transverse field Ising model: The transverse field 
Ising model (TFIM) is obtained by taking Si ■ Si = 3/4 in 
Eqn. (1), and setting J± = 0, H = Hx. In the vicinity of 
the TFIM's critical point (i.e. J z = H), this theory has a 
continuum representation as a free Majorana fermion: 20 

i r R v - - 

H = — dx-{i>d x i> + " * A V^- (4) 

27T Jo 2 

Here ip{x,t) and ip{x,t) are the right and left compo- 
nents of a Majorana Fermi field. The gap, A, of the 
fermions in the disordered regime (J z > H) is given 
by A - [J z - H). The Hilbert space of the the- 
ory on a periodic line of finite length R divides itself 
into two sectors: Neveu-Schwarz (NS) and Ramond (R). 
The NS-sector consists of a Fock space built with even 
numbers of half-integer fermionic modes, i.e. states 
of the form \p x ■ ■ -p 2 n)ns = apt ■ ■ ■a p2N \0) NS where a 
mode's momentum satisfies, pi — lixijii + 1/2)/ R, with 
rii an integer, while the R-sector consists of a Fock 
space composed of odd numbers of even integer fermionic 
modes, \k% ■ ■ ■ k 2 M+i)n = {a f kl ■ ■ ■ 4 2M+1 fc » = 

2irrii/R. The energy /momentum, E(pi)/P(pi), of a NS 
state, \pi ■ ■ -P2n)ns, is given simply by E(pj )/P{pj) = 
J2i e (Pi) I Pi where e(p) = y/p 2 + A 2 , with an 
identical relation holding for states in the R-sector. 

To compute the DSF of this model, we need access 
to the matrix elements, (l\S z \m), of the spin operator 
S z at finite R. These matrix elements, derived in Refs. 
(|20ll2ll ). only are non-zero when \l) and m) belong to 
different sectors. For such matrix elements we have 

R (k x ■ ■ ■ fc2M+i|^(0)|pi • • -P2n)ns = C R W_g{9 ki )g{0 Pj ) 

x n - ) n - e n ) n r 1 ^ - & 

i<j i<j i,j 

where p parameterizes the momentum pi via 
Pl = Asinh(0 Pi ), C R = 1 + 0(e- AR ), g{6) = 
(AJ?cosh(6»))- 1 / 2 , and f(6) = tanh(0/2). Note that up 
to exponentially small corrections, these matrix elements 



have the same functional form as at R = oo. The sole 
difference in the two cases is that at finite R, the mo- 
menta are quantized. This is a pattern that repeats itself 
for general integrable models, as emphasized in Ref. I22I 
and that we will exploit for our analysis of spin-1 chains. 

Crucially all of these matrix elements are finite, a con- 
sequence of working at finite R. The sole possible di- 
vergence comes from the term, Yii j / — 9pj), an d 
occurs as two momenta, ki and pj , approach one another. 
But as ki lies in the R-sector with integer quantization 
and pj lies in the NS-sector with half-integer quantiza- 
tion, the two are never exactly equal provided R is finite. 
In contrast, with R = oo the distinction between the R- 
and NS-sectors collapses (via a spontaneous Zi symme- 
try breaking). Concomitantly, S z has matrix elements 
where ki and pj may, in principle, be equal, and so which 
are infinite. By working at finite i?, we thus obtain a 
clean, unambiguous regulation of these infinities. 

Even though any given matrix clement is finite, we 
must still sum an infinite number of matrix elements in 
the Lehmann expansion of Eqn. §3§ in order to obtain 
the DSF, x 1 j a t finite T. To do so we employ a Dyson 
like equation by writing x 7 (w, q) in the form, 

X J (w, q) = D\uj, q)/(l - D t (oj, g)£>, q)). (6) 

D 1 ^, q) is the DSF in the absence of temperature in- 
duced interactions: D T (u),q) = 2e(q) / ((uj + iS) 2 + e 2 (q)) . 
As x^^il) has a well-defined low temperature expan- 
sion, so must S / (w, q): T, J (uj,q) — ^n^ni^'V)' where 
is at best 0(e~ n/3A ). We can readily compute £{. 
To do so we expand x 1 to first order in S{: X 1 — 
D 1 + (£> J ) 2 £{ +0(e- 2 P A ). We then compare this expan- 
sion to the expansion of x 1 i n terms of the Lehmann ex- 
pansion of Eqn. ©. To facilitate this we divide C I (x, r) 
into contributions coming from matrix elements with a 
fixed number of excitations on either side of the operator, 
S z , i.e. C ! {x,t) =Y,m,n C Iln( x > t ) where 

C m,n( x > t ) = ^2 \n(ki ■ ■ ■ k 2M +i\S z (0)\pi ■ ■ -P2n)ns\ 2 

k l'-' k 2M + l 

xe -0E(k j ) e -r(E(p i )-E(k j ))+ix(P( Pi )-P(k j )) _ / 7 N 

Kinematic constraints give that C l M N (u>,q) is of or- 
der e -/3max(2JVA-e(w)cj,(2M+i)A+0(-u;)w)_ Keeping terms 
to at least C(e _/3A ) and such that Imx 7 (— lj, q) = 
-Imx 7 (w,g), we reduce xH^^l) to X 1 = ^(Cox + C w + 
C\i + C21 + e-^^D 1 ). The final term, e^^D 1 , 
is a 'disconnected' contribution arising from C23 + C32, 
exactly cancelling off a similar contribution appearing 
in C21 + C\2- Comparing these two expansions for x 1 
gives us an expression for Si. First expanding out 
the partition function, Z 1 — J^^Lo^n where Zq = 1, 
Z[ = Y, P eR e ~ Mp) and generally Z^ is C( e -"' 3A ), and 
then noting that + C ! w = (1 - e^^ e{q) )D I , we obtain 
for £(, S{ = (C( 2 + C^KD 1 )- 2 (Z[ + e~^)(l - 
e -^w)(D 1 )- 1 . To then evaluate E{, we compute C[ 2 
and C21 numerically. As validation of our use of a finite 
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R regulation of the singularities in the matrix elements, 
S{ remains finite as R — » oo even though C[ 2 , C^ii and 
Z{ all diverge. 
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FIG. 1: Plots of the finite T DSF for the TFIM. 

In Fig. 1 we plot the resulting DSF, Imx 1 ^^ 
at a variety of temperatures. At T = 0.25A, the DSF 
is approximately Lorentzian, but as the temperature is 
increased to T = 0.6A, the lineshape develops a marked 
asymmetry. This asymmetry was also found to be present 
in a virial-like expansion of the finite T DSF. 23 We quan- 
tify the amount of asymmetry in the lineshape by com- 
puting the ratio, Alr, of the spectral weight to the left 
and the right of co = A. In the inset of Fig. 1 we compare 
our results for T = 0.5A (black solid curve) to those ar- 
rived at by using a semi-classical approach^ 3 - (blue curve) 
and to a Lorentzian approximation thereof (red dashed 
curve). We see that our computation yields asymmetries 
in the lineshape far in excess of those found in the semi- 
classical approach. While the semi-classics has a slight 
asymmetry at T — 0.5A, it is close to being Lorentzian. 

The origin of this discrepancy between the semi- 
classics and our treatment lies in two factors. The 
TFIM, as written in Eqn. ([4J, is relativistically invari- 
ant whereas the semi-classical model used in Refill only 
possesses Galilean invariance. It has already been noted 
that a relativistic dispersion relation, in comparison with 
a Galilean invariant one, better matches the measured 
lineshapeji^ A further difference arises in that the semi- 
classics for the TFIM is only strictly correct in the T — > 
limit. At finite T it misses corrections that here are en- 
coded in the form of the matrix elements (Eqn. ©). 
Spin-1 Heisenberg model: We now apply our ap- 
proach to the thermal broadening of the coherent mode 
in a gapped isotropic spin-1 chain (i.e. taking S • S = 2, 
J± = J z = J, and H = in Eqn. (ffj) . The isotropic 
spin 1-chain is given in the continuum limit by the 0(3) 
non-linear sigma models 



C = (2g) 1 (d x n ■ d x n + d T n ■ d T n) 



(8) 



The lattice spin operators, Si, are related to the con- 
tinuum fields by Si ~ (— l) l n(m ) + |n x d t n (with 



ao the lattice spacing). 2 - In this letter we will focus on 
the DSF near wavevector q = ir and so be interested 
in computing C(x,t) = (n z (x, r)n(O)). The spectrum 
and scattering matrix of the 0(3) nonlinear sigma model 
(NLSM) are known exactly. There are three elementary 
excitations, A\{9), a = x,y,z, forming a vector repre- 
sentation of 0(3). The excitations have a gap behav- 
ing as A ~ Je~ x / g . The excitations' energy and mo- 
mentum are parametrized in terms of the rapidity 9 via 
e(0) = A cosh(6>) and p{6) = A sinh(0). 

Like the TFIM, the eigenstates of the 0(3) NLSM can 
be delineated exhaustively in terms of multi-excitation 
states, i.e. \6\, a x ; • • • , 6 n , a n ) = A\ ± (6q) • • • A\ (0jv)|O). 
However, the matrix elements of these states involving 
the operator, n z , are considerably more complicated than 
those of the TFIM. But as we are working at low temper- 
atures, to compute the DSF, x i we W1 U only need re- 
course to matrix elements involving a maximum of three 
excitations (as with the TFIM). In infinite volume they 
are given by£>2£ 



(O|n°(O)|0,6) = 5 ab , 
{e u a 1 \n a (0)\e 3 ,a 3 - i e 2 ,a 2 ) 



x S aai 5 a 

where ip{0) 
0i-6 2 . 
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(9) 



9+in tanh 2 (fl/2) 
6{2iri+6) 9 ' 



9 — in, and 9\ 



As with the TFIM model, we work in finite volume. 
The sole effect of doing so upon the matrix elements, 
up to negligible e~ AR corrections, is to quantize the mo- 
mentum (i.e the 0's) with an attendant effect upon finite 
volume phase spaced 2 - Here however the quantization is 
more complex than that of the TFIM model. We must 
take into account the non-trivial interactions between the 
excitations and solve the Bethe ansatz equations. For the 
calculation at hand, we, at most, most solve the one and 
two particle Bethe equations for the states, |$i,ai) and 
|03) 03; 02, 02)- For the one particle state, we have the free 
quantization condition, 9\ = sinh _1 (27rn/i?) for some in- 
teger n. For the two particle case, 02/03's are quantized 
via e «A s inh(0 2 /9 3 ) _ e iS a {e 23 /e 32 )^ w ]j ere t ^ c n on-trivial 

phase, S a (9), marks the presence of interactions and de- 
pends upon the particular SU(2) representation (a = 
singlet/triplet/quintet) into which the two particle state 
falls. Because of the presence of interactions, the finite 
R matrix elements of the form (8±,ai |n a (0)|#3, a 3 ;9 2 ,a2) 
are never infinite as the 0's never coincide. We again see 
finite R provides a clean regulation of the singularities 
present at R = 00. 

The remainder of the calculation of x° 3 follows in 
exact analogy to the TFIM. x° 3 ( w , l) takes the form 
of Eqn. ^ with D° 3 = D 1 . In this case the ex- 
pansion of C° 3 (x,t) = {n z (x,T)n z (0)), C° 3 (x,t) = 



4 



^2m,n ^m,n( x ^ r )' appears as 

c£>(*,t) = l(0i,---,eMK(o)|ei;---;e^}| 2 

ei;- ; e M 

Xe -i3E(ej) e -r(E(e i )-E(e j ))+ix(p(e i )-p(e j ))^ ^ 

where 0„ = {9i,ai\. The partition function here admits 
the expansion, Z° 3 = 1 + Z? 3 + C(e- 2/3A ), and Z? 3 = 
3 E„ e/3Ac ° sh(e " ) With these redefinitions of C m ,n, D 
and Z n , takes the same functional form as Sf . 
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FIG. 2: Plots of the finite T DSF for the 0(3) NLSM. 

In Fig. 2 we plot the resulting finite T DSF for the 
0(3) NLSM. We again find that the lineshape is charac- 
terized by an asymmetry that grows with temperature. 
We observe that this asymmetry is far stronger than that 
of the semi-classical analysis (see Fig. 10 of Ref. [hi ) 
whose lineshape is well described by a Lorentzian even 
up to temperatures of T — 0.5A. The origin of this dis- 
crepancy between our treatment and the semi-classics is 



similar to that of the TFIM, i.e. finite energy effects 
in both the scattering of excitations and in the matrix 
elements involving the operator, n z . 

The asymmetry in the lineshape can be interpreted 
in terms of a temperature dependent gap, A(T)J^ The 
gap can be extracted as the location of the peak of a 
Lorentzian fitted to the asymmetric lineshape. We plot 
A(T) vs T in the left inset to Fig. 2 and compare our 
computations with the neutron scattering measurements 
performed on the spin chain, YBaNiO^, in Ref. (flU). 
We see good agreement. For other purposes, for example, 
the analysis of the three magnon scattering continuum 
in the spin chain compound CsNiCl^ for temperatures 
T > 0.25A, we suggest the fitting function, 

Imxi^q) = A/((u - e(q) - B) 2 + C) 1 "^- 1 ). (11) 

In the right inset to Fig. 2 we show that this function 
provides a good fit of our computed T = 0.4A lineshape. 

In conclusion we have presented a method by which 
the lineshape of the coherent mode of gapped spin chains 
can be determined at finite temperature. This method, 
employing a continuum integrable representation of the 
chains, works with finite length, R, systems so as to 
circumvent infinities in matrix elements that appear at 
R = oo. It further employs a Dyson-like resummation 
of the matrix elements appearing in a Lehmann expan- 
sion of the DSF. The primary conclusion drawn from this 
analysis is that the lineshape of the mode is asymmetric 
with an asymmetry increasing with temperature. 
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